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The reliability of Langmuir probe measurements for plasma-turbulence investigations is studied 
on GEMR gyro-fluid simulations and compared with results from conditionally sampled I-V char- 
acteristics as well as self-emitting probe measurements in the near scrape-off layer of the tokamak 
ASDEX Upgrade. In this region, simulation and experiment consistently show coherent in-phase 
, fluctuations in density, plasma potential and also in electron temperature. Ion-saturation current 

measurements turn out to reproduce density fluctuations quite well. Fluctuations in the floating 
potential, however, are strongly influenced by temperature fluctuations and, hence, are strongly 
distorted compared to the actual plasma potential. These results suggest that interpreting floating 
■ as plasma-potential fluctuations while disregarding temperature effects is not justified near the sep- 

aratrix of hot fusion plasmas. Here, floating potential measurements lead to corrupted results on 
the ExB dynamics of turbulent structures in the context of, e.g., turbulent particle and momentum 
transport or instability identification on the basis of density-potential phase relations. 
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I. INTRODUCTION 



In the boundary region of magnetically confined fusion plasmas, micro instabilities determine a major part of 
transport losses across the separatrix into the scrape-off layer (SOL) and, therefore, play a key role for the global 

■ confinement quality. In order to understand the mechanisms degrading or even improving confinement via vortex-flow 
interactions, the involved dynamical processes are addressed by measurements of, e.g., correlations between density 

■ (n) and plasma-potential (^pi) fluctuations as well as average radial turbulent transport of particles (r = (nv r ) 
with radial ExB velocity fluctuations v r — Eg/B) and momentum (perpendicular Reynolds stress R s = (vgv r )), see 
e.g. [TI-fToj among others. 

The acquisition of local plasma-potential fluctuations as required for these studies is mostly based on Langmuir 
probe measurements. The measured floating potential of the probe is related to the plasma potential via 



i fl = $ pl - a— (i) 

On 

with A«3. Electron-temperature fluctuations (T e ) are rarely available in experiments and usually neglected such that, 
for simplicity, $a = <J> p i is used. The simplifying assumption that temperature fluctuations are negligible has been 
proven valid in a toroidally confined low-temperature plasma |Tlj], where temperature gradients were much smaller 
than the density gradient. However, in the boundary of several high-temperature plasmas, significant temperature 
fluctuations have been found by fast swee ping Langmuir probes [1214151 1, tripple probes [lj, EE E3]i swept double 
probes EH , a harmonic probe technique 0, Il9( and Thomson scattering [J, [20| • 
. £^ | In this work, the influence of temperature fluctuations on Langmuir probe measurements of density and potential 
fluctuations in a high-temperature plasma is investigated. To this end, simulations of plasma-edge turbulence are 
carried out with the gyro-fluid code GEMR. Synthetic probe data from the simulations are analyzed with respect 
to temperature effects. The results are compared to probe measurements near the last closed flux surface (LCFS) 
in the tokamak ASDEX Upgrade [Uj]. For the comparison, two independent sophisticated approaches are employed 
to resolve T e experimentally: one is based on self-emitting probe measurements, the other on conditionally sampled 
probe characteristics. It will turn out that temperature fluctuations near the LCFS of high-temperature plasmas 
strongly distort floating-potential measurements, and an interpretation of $a as $ p i is tenuous. 

This paper is organized as follows: A short description of the GEMR code is given in Sec. HH Section Hill presents 
the analyzes of the synthetic probe data from the simulations and temperature effects will be identified. The re- 
sults are compared with experimental data from both, self-emitting probe measurements and conditionally sampled 
probe characteristics, in Sees. [TV] and E] respectively. In Sec. |VH the results are discussed with respect to possible 
consequences for the interpretation of probe data and a conclusion is given. 
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II. GEMR SIMULATIONS 



The three-dimensional gyrofluid model GEMR [22| solves the first six moments of the gyrokinetic equation (density, 
parallel velocity, parallel and perpendicular temperature and associated parallel/parallel and perpendicular/parallel 
heat flux) for ions and electrons. The corresponding conservation equations are derived using a consistent treatment of 



the energy conservation 23] . The plasma species are connected through the quasi-neutrality condition obtained from 
the gyrokinetic polarization equation [24j and the induction obtained from Ampere's law. The number of moments 
retained in the model enables GEMR to capture typical tokamak edge regimes, where realistic density and temperature 
gradients yield a mixture of drift-Alfven and ion-temperature-gradient (ITG) turbulence. Collision frequencies may 
reach values below the Braginskii collisional regime. The gyrofluid formulation guaranties a proper presentation of 
the unit order finite Larmor radius effects in the drift wave component. 

GEMR is a global and not a flux-tube code. It allows the geometry to vary both radially and poloidally. A circular 
cross-section with toroidal axisymmetry is assumed. The coordinate system is aligned with the equilibrium magnetic 
field to take computational advantage of the strong spatial anisotropy of magnetised plasmas. Global consistency in 
the angles is ensured, even for toroidally truncated domains 25] . The present simulations were performed in a quarter 
of the torus with a resolution of 128 x 512 x 16 in the x, y and s directions, respectively. The coordinate system 
consists of the flux label x, the field line label y and the coordinate parallel to the magnetic field s. It is linked to the 
polar coordinates (r, 9, (f) by the unit Jacobian transformation rules 

x = 2ir 2 R r 2 y k = (q s 6 - - a k ) / (2tt) s = 9/ (2tt) (2) 

with a toroidal shift a k (r) yield by the shifted metric procedure due to the shearing of the magnetic field 26] . ASDEX 
Upgrade parameters were chosen for the tokamak major and minor radii of Rq = 1.65 m and tq = 0.5 m, respectively. 
The safety factor was set to q s — 4.6 and the magnetic field to B — 2.6 T. The inner half of the radial box corresponds 
to a 3 cm wide edge region and the outer half corresponds to the SOL, with a limiter cut located at the bottom of the 
flux surfaces. Within the local approximation of homogeneous normalization parameters, the code gives a consistent 
description of the turbulence from the confinement region across the LCFS into the SOL. In the SOL, the field lines 
intersect material plates, allowing the existence of modes with parallel and perpendicular wave numbers k\\ = and 
k± =/= 0, respectively. This changes substantially the behavior of heat and particle fluxes down the gradients [27|. The 
SOL is implemented with a change in the parallel boundary conditions according to a linearised Debye sheath model 
[13, HH, consistent with the formulation of GEMR, which retains only quadratic nonlinearities. 

The zonal profiles of the state variables (densities, temperatures and field potentials) are evolved together with the 
fluctuations. This accounts for the vigor of SOL turbulence. The linear drive terms are included in the gradients 
of the evolving profiles. Profile maintenance is achieved with source/sink zones at each radial boundary. The zonal 
components of the densities and temperatures are feedback dissipated towards an initially specified profile [28]. The 
constant density Hq and temperature T used for normalization in the simulations reflect typical ASDEX Upgrade 
L-mode parameters, namely To = T, = T e = 60 eV and no — rii — n e = 1.2 x 10 19 m~ 3 . The radial decay lengths 
of temperature (Lt) and density (L n ) have been chosen as Lt = L n /2 — 3 cm. The deuterium ion to electron mass 
ratio was rrii/m e = 3670 and the effective ion charge was Z e s = 2. The following section compares 4 ms of turbulent 
fluctuations from a saturated phase of the GEMR simulation with simultaneous Langmuir multi-probe measurements 
at the LCFS in the outboard midplane. 

III. SYNTHETIC-PROBE RESULTS 

Synthetic Langmuir probes have been implemented in GEMR to extract measurements from the simulated data. 
The ion-saturation current fluctuations are calculated according to the linearized equation 

T ~A lfe \ 

hat ~ A p en c s0 — — (3) 

y n 2 T J 

assuming T e = T,;, to be consistent with the GEMR model. Electron density (n e ) and temperature fluctuations (T e ) 
evolve with time, while uq and To are constant normalization parameters. The current is averaged over the parallel 
projection area (A p ) of an ASDEX Upgrade pin probe. The constant ion-sound velocity c s q is calculated using To. 
The results do not change significantly if the common expression 



0.6lA p enJ^^ (4) 
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FIG. 1: Comparison of simulated plasma fluctuations with synthetic probe data. Shown are details of the raw time series (left) 
and the conditionally averaged fluctuations (right) at the LCFS. Top: density and ion-saturation current (dashed) fluctuations; 
Middle: plasma and floating (dashed) potential fluctuations; Bottom: electron temperature fluctuations from GEMR. 



is used instead. The floating potential fluctuations are approximated by Eq. ([!]) with A = 3.4 for a cold Langmuir 
probe in a deuterium plasma. Different models predict similar values for A 29] . The differences are small, especially 
when compared with the experimental uncertainties concerning secondary electron emission, ion temperature and 
effective particle collection areas [30]. However, the exact value is not decisive for the nature of turbulent fluctuations, 
even if the absolute values change. 

The left-hand side of Fig. [1] compares the simulated plasma fluctuations n e (top), <E> p i (middle) and T e (bottom) in a 
100 /is time window with the synthetic probe data of 7 sa t and <&fl. All signals are taken from the LCFS at the outboard 
midplane of the circular simulation domain. Good agreement is found between fluctuations of plasma density and 
ion-saturation current. This indicates a minor influence of the temperature fluctuations and justifies the common 
experimental assumption J sat oc n e in the simulation. In contrast, no agreement is found between fluctuations of 
plasma and floating potential (Fig. [I]:). $a rather follows the inverse temperature behavior as shown in Fig.QJ, which 
according to Eq. (JT|) indicates the dominance of T e in the measured floating potential. 

In order to study the correlation between the fluctuating quantities systematically, the coherent part of the fluc- 
tuations is extracted by means of the conditional averaging technique [3l|. The ion-saturation current was used as 
reference signal with a trigger threshold of 2a (upper horizontal dashed line in Fig. [TJl) , where a is the standard 
deviation. The right-hand side of Fig. [1] shows the conditional average of 660 independent series of fluctuations in 
density and ion-saturation current (top), plasma and floating potential (middle) and temperature (bottom). The 
average temporal behavior of coherent structures in the synthetic ion-saturation current and in the density is again 
found to agree very well. However, the synthetic floating potential is clearly anti-correlated with the simulated plasma 
potential fluctuations (Fig.QJl). This can be attributed to significant coherent temperature structures as depicted in 
Fig. [Tf, which show up in phase with structures in n and $ p i. 

The analyses carried out on data from turbulence simulations suggest that close to the LCFS in high-temperature 
plasmas, floating-potential fluctuations measured with probes are strongly affected by temperature fluctuations and, 
therefore, do not even qualitatively reflect the plasma potential. This reflects the experimental situation if two 
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FIG. 2: Photo (left) and sketch (right) of the multi-pin probe head as seen from the plasma in the tokamak ASDEX Upgrade. 
The shaded areas in the sketch are retracted by 4 mm and 8 mm (probe 13). 

assumptions hold: (I) The turbulent fluctuations calculated by the gyro-fluid simulation are comparable to plasma 
fluctuations in the boundary of a real tokamak. (II) In hot magnetized plasmas, the behavior of Langmuir probes 
is sufficiently well described by the Debye sheath model. Hence, the presence of significant temperature fluctuations 
would have severe consequences for the interpretation of fluctuations from probe measurements as will be discussed in 
Sec. I VII In the next two sections, it will be shown that the above assumptions are indeed supported by experimental 
observations. 



IV. EMISSIVE-PROBE MEASUREMENTS 



In order to test the validity of the simulation results from the previous section, probe measurements were carried 
out in the tokamak ASDEX Upgrade [2ljj . The experiment was conducted in an L-mode deuterium plasma, kept 
in lower single null divertor configuration with a constant toroidal magnetic field of 2.5 T (clockwise) and a plasma 
current of 800 kA (counter clockwise from top view) . The auxiliary electron cyclotron heating power was 600 kW and 
the edge plasma density 4.0 x 10 19 m~ 3 . A Langmuir probe array penetrated the plasma horizontally 0.31m above 
the outer torus midplane [32j . 8 free standing carbon pins on the top level measured alternatively floating potential 
and ion-saturation current, as shown in Fig. [2] (white background). The I sa t probes were biased with —180 V and 
the floating-potential signals were amplified directly at the manipulator exit to minimize low-pass filtering by cable 
capacity and plasma impedance. The data were sampled at 2 MHz with 14 bit resolution. For the investigation of 
only the fast dynamics upon the stationary background, fluctuations below 10 kHz have been digitally removed from 
the signals. 

For comparison with simulation results from the previous section, fluctuations in the floating and plasma potential, 
density and temperature are deduced from the probe signals in two different, independent ways: The first method 
- presented in this section - involves direct measurements of the plasma potential with an emissive probe. The 
next section is dedicated to the second method, which is based on conditionally sampled probe characteristics. Both 
methods are applied to data from the same discharge. Figure [3] shows raw data from three Langmuir probes during 
the second stroke of discharge #26530. The ion-saturation current (J sa t, top) increases and the floating potentials 
($fl, bottom) decrease while the probes approach the LCFS (t < 3.57s). At t » 3.582s, probe number 10 (red) 
reaches a critical temperature due to the heat flux from the plasma onto the probe and starts to emit electrons. The 
emitted electron current is space charge limited after a fast transition within 2 /is. Now. the self-emitting probe floats 
at the potential 

$ sc = $ P i - 0.6—, (5) 

e 

which corresponds to the sheath edge separating sheath from presheath. This was shown in previous investigations 
with self-emitting probes in ASDEX Upgrade J33] ■ Probe number 9 is less heated and measures the floating potential 
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FIG. 3: Measurements of ion-saturation current (top) and potential (bottom) during a stroke of Langmuir probes in ASDEX 
Upgrade. Due to the high heat load in the plasma, one probe emits electrons during the outward motion and measures a more 
positive potential (red line). 



of a cold probe 

$ fl = * pl - 3.4^ (6) 

e 

during the entire discharge. 

Combining Eqs. ([5]) and © the electron-temperature dynamics can be deduced from simultaneous measurements 
of $fl and <I> se with a temporal resolution of 0.5 fj,s. Furthermore, using the result for the temperature and assuming 
T e p» Ti, the density can be calculated according to Eq. ^ from nearby measurements of the ion-saturation current. 
Similar to Fig. [TJ Fig. [4] shows short fluctuation time traces of the results on the left-hand side and the conditional 
average of 28 independent large amplitude events with J sat > 2a on the right-hand side. The stationary fluctuations 
are taken from a 4.1ms time interval at t = 3.625 s, which corresponds to the radial range 14 — 18 mm outside the 
LCFS. An anti-correlation is observed between floating and plasma potential (Fig. 0J1), similar to the simulation 
result in Fig. [TJi. Electron temperature fluctuations with a significant amplitude (Figs. [4^, f) are found in phase 
with plasma-potential fluctuations. Fluctuations in the density are almost coincident with those in J sat (Figs. 
b), i.e. the temperature fluctuations do not carry significant weight. In comparison, a good agreement of simulation 
and experimental results is found, which confirms the crucial role of temperature fluctuations in interpreting probe 
data. It can be shown, that the experimental findings are not distorted by the probe separation, which introduces 
additional phase delays between the signals due to the perpendicular propagation of turbulent structures. Structures 
propagating with typically 1 km/s would show a phase shift of 5 [is over a distance of 5 mm. This might explain a small 
phase delay, but does barely account for the anti-correlation of the floating and the plasma potential. Nevertheless, a 
different method, which does not suffer from the probe distances, is used in the next section to verify the experimental 
results given here. 



V. CONDITIONAL SAMPLING OF PROBE CHARACTERISTICS 



In this section, the previous results from emissive probes are verified via an alternative method, which resolves 
fluctuations in potential, electron temperature and density quasi-instantancously at the same position. This method 
has be en p roposed recently and is based on the conditional sampling of characteristics from a slowly swept Langmuir 
probe [34(. Here, it is applied for the first time to the high-temperature plasma in ASDEX Upgrade. Figure [S] shows 
time windows of an ion-saturation current measurement (7 sat , top) together with the current (/, middle) and voltage 
(V, bottom) of a slowly swept Langmuir probe nearby the / sa t measurement. Complete probe characteristics are 
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FIG. 4: Comparison of measured plasma fluctuations from cold and hot Langmuir probes in the SOL of ASDEX Upgrade. 
Detail of time series (left) and conditionally averaged fluctuations (right). Top: density and ion-saturation current (dashed); 
Middle: plasma and floating (dashed) potentials; Bottom: electron temperature from Eqs. ([5]) and ((6]). 
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FIG. 6: Results from fitting conditionally sampled I-V characteristics close to the LCFS of AUG. Left: density and ion- 
saturation current (dashed); Middle: plasma and floating potential (dashed); Right: electron temperature from fits to probe 
characteristics. 



swept with a rate of 1 kHz. After conditional sampling, a full characteristic is available every microsecond within 
a short time window around the trigger condition of 1.5a amplitude in the J sat signal of probe number 2. Around 
each trigger event at ti, a short I-V trace is collected from the neighboring probe number 6. From several trigger 
events, a complete I-V characteristic can be reconstructed at each time lag At with respect to the trigger times U, 
since the trigger is not phase locked to the sweeping bias voltage. In the present case, 510 trigger events have been 
detected within a 40 ms time interval. During this time, the probe was located in the SOL 8 — 15 mm outside the 
LCFS. The temporal resolution is reduced from 0.5 to 1 /is, i.e. 1020 I-V pairs are combined to one conditionally 
sampled I-V characteristic. This is sufficient to obtain reliable fits to the characteristics providing T e , <frfl and I sa t 
every microsecond. 

Figure [5] shows the fluctuations obtained from fits to the conditionally sampled I-V characteristics, together with 
plasma density and potential calculated according to the Eqs. Q (T e « Ti) and ©, respectively. The results 
are qualitatively the same as those from emissive-probe measurements: Density, plasma-potential and temperature 
fluctuations are essentially in phase. The ion-saturation current does not substantially deviate from the density. 
The floating potential, however, is seriously affected by the temperature fluctuations resulting in an anti-correlation 
between floating and plasma potential. Moreover, the experimental results confirm once more those from the numerical 
simulations in section UTTl This means that (I) the plasma fluctuations in the boundary region of ASDEX Upgrade are 
sufficiently well represented by the gyro-fluid simulation and (II) the Debye sheath model is appropriate to describe 
Langmuir probes under the present conditions. It turns out clearly that great care has to be taken in interpreting 
data from Langmuir probes when strong temperature fluctuations are present. 



VI. DISCUSSION AND CONCLUSION 



The objective of this paper was to study the influence of electron temperature fluctuations on turbulence investi- 
gations with Langmuir probes in high-temperature plasmas. To this end, gyro-fluid simulations of plasma turbulence 
were carried out for fusion edge plasmas comparable to that in ASDEX Upgrade. Synthetic Langmuir probes were 
implemented to mimic the experimental diagnostics. From the simulations, coherent structures in the fluctuations 
of the electron temperature could be expected in the SOL close to the LCFS of a real tokamak. For comparison, 
Langmuir-probe measurements were carried out in the edge region of ASDEX Upgrade. The simulation results could 
be confirmed in two independent experimental approaches, one based on self-emitting probe measurements and the 
other on conditionally sampled probe characteristics. 

In investigations on plasma dynamics, probe measurements of the ion-saturation current and the floating-potential 
fluctuations are widely used interchangeably for the density and the plasma-potential fluctuations, respectively. In this 
work, it was shown that coherent temperature fluctuations seriously alter the potential measurements. In simulation 
and experiment it was found that close to the LCFS, coherent density, plasma-potential and electron-temperature 
fluctuations are essentially in phase. Density fluctuations are reasonably well represented by the ion-saturation current, 
i.e. the influence of T e is marginal. Plasma and floating potential, however, show a completely different behavior: 
they are found to be anti-correlated as a consequence of the significant T e contribution. An influence of T e on <I>fl is 
also observed in a numerical investigation [35[ . 
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This has far-reaching consequences for probe based turbulence investigations. The phase relation between density 
and plasma potential is an indicator for the turbulence-driving instability. It determines whether a density perturba- 
tion is stable or unstable depending on the ExB advection associated with the corresponding potential perturbation. 
Hence, the cross phase determines the level and even the direction of turbulent particle transport. This work shows, 
that the phase between plasma and floating potential is governed by the potential-temperature cross phase. In turn, 
the floating potential does not reflect the plasma potential. Deduced quantities like density-potential cross phases 
and turbulent transport as well as Reynolds stress involving electric field estimates would be distorted if based on 
floating potential measurements. 

For example, the floating potential in Fig. [6p suggests that the density is accompanied by a negative monopole 
structure in the potential. In fact, the plasma potential shows a positive monopole structure, as expected from the 
high parallel electron mobility. As a result the ExB drifts are inverted and indicate an opposite vorticity. From 
floating-potential measurements, the dynamics of coherent potential structures would hence be misinterpreted. In 
particular, the average radial turbulent transport F — (nv r ) with v r = Eg/B and Eg « —(^a,2 — $a,i)/A (A is the 
poloidal probe distance) would point into the wrong direction. Similarly, estimates of the perpendicular Reynolds 
stress (R B = (vgv r )) as a measure of vortex tilting can be affected qualitatively if temperature fluctuations deform the 
topology of the floating compared to the plasma potential. From the formal construction of R s which is quadratic in 
the electric field, it can be seen, however, that the polarization alone is not relevant. I.e., R s measurements near the 
LCFS might be at least qualitatively valid. Moreover, it is worth noting that in the present study, the fluctuations in 
the SOL near the LCFS of hot fusion plasmas reveal drift-wave like features, i.e., a density-potential cross phase close 
to zero. In this region, density blobs come along with a monopole potential rather than a delayed dipole structure, 
which would be associated with an interchange mechanism behind the blob dynamics. The results are consistent 
with previous plasma-potential measurements close to the confined plasma of TJ-K (3(|, DIII-D Q and ASDEX 
Upgrade [37| . The coincidence of density and electron-temperature fluctuations reflects the transport of heat into the 
SOL and is consistent with findings from TEXT, Phaedrus-T, TJ-I, W7-AS and Repute-I [38| . 

In conclusion, plasma fluctuations in the SOL near the LCFS of hot plasmas in ASDEX Upgrade are well described 
by GEMR gyro-fluid simulations. In this region, measurements and simulation showed the existence of coincident, 
coherent structures in density, potential as well as electron-temperature. The temperature fluctuations are large and 
coherent enough to significantly distort floating-potential measurements with Langmuir probes for interpretation as 
plasma-potential fluctuations. At large temperature amplitudes, the distortion of the floating potential generally 
depends on the phase and coherence of the temperature-potential coupling. Hence, great care must be taken when 
interpreting floating-potential measurements with respect to, e.g., cross-phase relations and turbulent transport of 
particles or momentum. Results obtained in the SOL near the LCFS of hot fusion plasmas similar to ASDEX Upgrade 
are likely to be corrupted if the effect of temperature fluctuations is disregarded. In this work, methods based on self- 
emitting probe measurements or conditionally sampled probe characteristics proved useful to test for such temperature 
effects. 
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